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This study investigates the relative influence of biotic and abiotic factors on community 
dynamics using an integrated approach and highlights the influence of space on genotypic 
and phenotypic traits in plant community structure. We examined the relative influence 
of topography, environment, spatial distance, and intra- and interspecific interactions on 
spatial distribution and performance of Boechera stricta (rockcress), a close perennial 
relative of model plant Arabidopsis. First, using Bayesian kriging, we mapped the 
topography and environmental gradients and explored the spatial distribution of naturally 
occurring rockcress plants and two neighbors, Taraxacum officinale (dandelion) and 
Solidago missouriensis (goldenrod) found in close proximity within a typical diverse 
meadow community across topographic and environmental gradients. We then evaluated 
direct and indirect relationships among variables using Mantel path analysis and developed 
a network displaying abiotic and biotic interactions in this community. We found 
significant spatial autocorrelation among rockcress individuals, either because of common 
microhabitats as displayed by high density of individuals at lower elevation and high 
soil moisture area, or limited dispersal as shown by significant spatial autocorrelation of 
naturally occurring inbred lines, or a combination of both. Goldenrod and dandelion density 
around rockcress does not show any direct relationship with rockcress fecundity, possibly 
due to spatial segregation of resources. However, dandelion density around rockcress 
shows an indirect negative influence on rockcress fecundity via herbivory, indicating 
interspecific competition. Overall, we suggest that common microhabitat preference 
and limited dispersal are the main drivers for spatial distribution. However, intra-specific 
interactions and insect herbivory are the main drivers of rockcress performance in the 
meadow community. 



Keywords: Bayesian kriging, competition, correlogram. Mantel tests, path analysis, spatial interaction, spatial 
pattern 



INTRODUCTION 

Spatial patterns are a crucial, but often overlooked, component to 
understanding the factors and processes that structure plant com- 
munities (Levin, 1992; Jeltsch et al., 1999; Mclntire and Fajardo, 
2009). Spatial distribution and performance (i.e., growth and 
reproduction) of any plant species depends on the ability to 
cope with environment, most notably topography, soil properties, 
and moisture availability (Goslee et al., 2005), intraspecific geno- 
type variation of plant species (Crutsinger et al., 2006; Lankau 
and Strauss, 2007), intra- and interspecific plant-plant interac- 
tions (Callaway and Walker, 1997; Holzapfel and Mahall, 1999; 
Pugnaire and Luque, 2001; Brooker et al., 2008; Genung et al., 
2012), and insect herbivory (Marquis, 1992; Bloom et al., 2003; 
Becerra, 2007). Investigating the relative importance of these 
controlling factors is critical to understanding the underlying 
processes that structure plant communities. This study describes 
the first application of an integrated approach using molecular, 
ecological, and statistical tools to quantify the relative influence 
of biotic and abiotic factors on the spatial distribution and 



performance of Boechera stricta (rockcress), an emerging ecolog- 
ical model plant species (Rushworth et al., 2011). Such species 
show great promise for understanding genetic controls on eco- 
logically important traits (Song and Mitchell-Olds, 2011) and 
provide opportunity to explore underlying mechanisms in nat- 
ural populations, compared to artificial experimental settings. 
Rockcress is widely distributed in the western United States (Song 
and Mitchell-Olds, 2011) and shows local adaptations in diverse 
ecological habitats (Mitchell-Olds, 2001; Knight et al, 2006; Song 
et al, 2006), making it an ideal species to study the biotic and 
abiotic controls on species distribution and growth in natural 
settings. 

Combinations of complementary mathematical and statisti- 
cal techniques have been used in recent studies to investigate 
links between observed spatial patterns and underlying ecolog- 
ical processes [e.g., process models and geostatistics (Loranty 
et al, 2008); multiple ordination and geostatistics (Wagner, 
2003); Structural Equation Modeling and Bayesian statistics 
(Arhonditsis et al., 2006); mixed models and Bayesian kriging 
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(Smithwick et al., 2012); and Mantel tests and path analysis 
(Goslee et al., 2005)]. The integration of these tools can pro- 
vide information about spatial structure of variables and potential 
underlying ecological processes. Here, we integrated Bayesian 
kriging (Diggle et al, 1998; Diggle and Ribeiro, 2002), a novel 
data-model fusion approach for spatial interpolation of topog- 
raphy and environmental data across the local landscape, and 
Mantel path analysis (Mantel, 1967; Leduc et al., 1992; Goslee 
et al, 2005) to infer the relative importance of microhabi- 
tat preference, limited dispersal, competition, and herbivory in 
determining the spatial distribution and performance on a local 
scale. 

First, we examined the spatial distribution of rockcress indi- 
viduals and the density of inter-specific neighbors across the 
topography of the local landscape, which included gradients 
of soil moisture, vapor pressure deficit, and topology (mea- 
sured as elevation). Second, we determined the spatial struc- 
ture of environmental and intra- and interspecific variables 
(see Methods for detailed description of these variables) using 
piecewise Mantel correlograms (Goslee and Urban, 2007). The 
Mantel correlogram is particularly useful for studying eco- 
logical patterns in count data because it provides the spatial 
information in terms of distance apart rather than geograph- 
ical location using dissimilarity-based analysis (Urban et al., 

2002) . Third, we evaluated plausible hypotheses on underly- 
ing processes using Mantel path analysis in the presence and 
absence of space (Goslee et al., 2005). We asked: (1) what gov- 
erns the spatial distribution and performance of rockcress on 
a local scale? and (2) what is the relative importance of fac- 
tors governing the spatial distribution and performance at this 
scale? 

MATERIALS AND METHODS 
STUDY SITE 

We studied spatial distribution and performance of rockcress 
[Boechera stricta (previously Arabis drummondii (Al-Shehbaz, 

2003) ] plants occurring within a conspicuous highly populated 
local area in the northern Black Hills, South Dakota, USA, 
(Lat: 44.403611, Long: -103.938403, elevation 1365 m) during 
the summer of 2004. Rockcress plants grew within a 40 x 50 m 
meadow area surrounded by ponderosa pine (Pinus ponderosa), 
aspen (Populus tremuloides), birch (Betula spp.), and burr oak 
(Quercus macrocarpa). Other common plant species that inhab- 
ited the meadow were goldenrod (Solidago missouriensis), dan- 
delion (Taraxacum officinale), mouse-ear chickweed (Cerastium 
vulgatum), snowberry (Symphoricarpos albus), and sedge (Carex 
spp.). The relief between the highest and the lowest measured 
point across the study area was 9.55 m. 

ENVIRONMENTAL DATA 

A hand-held probe (Theta Probe, Delta-T Devices) was used 
to measure surface (0-6 cm) soil moisture and a hand-held 
sling psychrometer (Bacharach 12-7012; Bacharach, Inc., New 
Kensington, PA) was used for vapor pressure deficit. These mea- 
surements were taken at 71 randomly selected locations once to 
represent the spatial gradient of environmental variables across 
the meadow. 



SAMPLING 

Every rockcress individual (n = 234) in the study area was geolo- 
cated and sampled for growth (diameter of basal rosette, rosette 
diameter) and reproductive fitness (number of reproductive 
stalks, stalk number; reproductive stalk diameter, stalk diame- 
ter; and number of fruits, fruit number). In addition, amount 
of herbivore damage (percent area consumed) on basal rosette 
leaves [rosette herbivory) and reproductive stalk leaves {stalk her- 
bivory) were recorded. Herbivory was quantified for each plant 
by counting the number of holes (>1 mm) that the flea beetles 
chewed on the leaves of basal rosette and reproductive stalks. 
To determine the influence of neighboring species on rockcress 
growth and performance two plant species, goldenrod and dan- 
delion were selected due to their abundance in the landscape in 
close proximity to rockcress as compared to other species. The 
number of goldenrod (goldenrod density) and dandelion (dan- 
delion density) plants within a 15 cm radius of each rockcress 
plant was recorded. The 15 cm radius was an appropriate dis- 
tance because of the small size and high density of the meadow 
plants. 

Additionally, to distinguish between environmental and 
genetic causes for spatial patterns, we genotyped 142 randomly 
selected rockcress plants using seven polymorphic microsatellite 
loci and STRUCTURE software (Pritchard et al, 2000), avail- 
able at http://pritchardlab.stanford.edu/structure.html. Briefly, 
this technique identifies groups of relatives, which in this case 
are self-fertilizing (inbred) lineages (hereafter, line), or in other 
words, overlapping-generation dependents from common ances- 
tors. Please see detailed description of this technique in Siemens 
et al. (2009). 

STATISTICAL PROCEDURES 
Bayesian kriging 

Maps of environmental variables were generated using spatial 
interpolation in a Bayesian framework. A fully probabilistic 
Gaussian spatial process model (Diggle et al., 1998; Diggle and 
Ribeiro, 2002) was used for Bayesian kriging, which assumes that 
conditional on a Gaussian underlying process, S the observed 
data, Yi: i = 1, . . . n are independent with a distribution in the 
exponential family. A brief summary of the modeling approach 
is given below; detailed information can be found in Diggle and 
Ribeiro (2002). The model can be expressed in a hierarchical 
framework as follows: 

Level 1 : Yi\S ~ N(B(xi) + S(xi), r 2 ) 

Level 2 : S(xi) ~ N(0, o 2 R(h; O)) (1) 

Level 3 : prior( / S, a 1 , O, r 2 ) 

The first level describes a spatial linear trend (B = trend param- 
eter) based on spatially referenced explanatory variables, r 2 
(nugget) represents measurement variability and/or spatial vari- 
ation below the sampling grain. The second level describes a 
stationary Gaussian spatial process (S(xi)) with mean = 0, vari- 
ance = a 2 and correlation function R(h; <J>), where <J> is cor- 
relation parameter (range of spatial autocorrelation = 3<J>) and 
h is lag distance (vector distance between two locations), and 
the third level specifies the prior for the model parameters. We 
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chose an exponential correlation function to quantify spatial 
autcorrelation: 

R(/i; <5) = exp(- h/<&) (2) 

The mean and variance of topography and environmental data 
were estimated at individual locations from the predictive dis- 
tribution using the krige.bayes function of geoR library (Ribeiro 
and Diggle, 2001) in R version 2.15.2 (R Development Core 
Team, 2012). This algorithm samples a parameter value from a 
discrete posterior distribution computed from joint distribution 
of parameters and priors. We assumed a constant trend mean 



model and chose a sensible interval of values for each param- 
eter considering the study site to generate a multidimensional 
parameter [<J>, a 2 , and r 2 .rel (relative nugget = r 2 /cr 2 )] grid. 
Please refer to Appendix for exact intervals for individual param- 
eters. Flat priors (see Figures A1A-C for an example of prior 
and posterior distributions) were chosen for <t>, and r 2 . rel, and 
a reciprocal prior for a 2 . The sampled parameter value is then 
attached to [ft |Y, <t>, a 2 , r 2 .rel] and a realization is obtained 
from the predictive distribution at an unsampled location. A large 
sample size was generated by repeating this process several times 
which allowed a stable estimation of the underlying distribution. 
The mean and the variance of the predictive distribution were 
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FIGURE 1 | Map of spatial distribution of (A) Boechera stricta plants, 
each circle indicates the presence of one plant and diameter of the 
circles corresponds to the rosette diameter classes from 1-24, 25-45, 
and 46-60 cm. Rosette diameter zero is marked with an x indicating 
previously present plant which is currently dead; (B) a random sample (142 



individuals) of B. sticta plants showing 6 putative naturally occurring 
near-inbred lines (line 1-6) (CD) number of Solidago and Taraxacum plants 
within a 15cm radius around each B. stricta plant, diameter of the circles 
correspond to the number classes 1-20, 21-40, and 41-60. (E,F) Kriged maps 
of elevation and soil moisture for the study site. 



www.frontiersin.org 



July 2014 | Volume 5 | Article 348 | 3 



Naithani et al. 



Spatial distribution and performance of Boechera striata 



computed at unsampled spatial locations using 100,000 posterior 
draws. Leave-One-Out cross-validation (Figure Al in Appendix) 
was used for model validation. The maps of environmental vari- 
ables were used to obtain data at individual plant's location to 
conduct the Mantel correlation analysis. 

The piecewise Mantel correlogram and Mantel test 

A correlogram in this context is a plot of the spatial autocorre- 
lation with lag distance, and the Mantel test can be used to test 
the significance of the correlations. A simple Mantel test (Mantel, 
1967; Goslee and Urban, 2007) was performed using Euclidean 
distance to explore variation in spatial structure of topogra- 
phy (local vertical relief measured as elevation), environmental 
drivers (soil moisture and vapor pressure deficit), fitness param- 
eters [growth (rosette diameter) and reproduction (stalk numbers, 
stalk diameter, fruit number), and line] of rockcress, density of 
neighboring plants (goldenrod density and dandelion density), 
and herbivory on rockcress leaves (rosette herbivory and stalk 
herbivory). Data were standardized before calculating Euclidean 
distances (as recommended by Goslee, 2010) and the error was 
computed using 10,000 permutations. Autocorrelation range was 
determined by looking at the significance value (P = 0.05) of the 
individual bins of lag distance at every 2.5 m in the correlograms. 
A pairwise correlation between all the variables was determined 
using simple Mantel test that guided the path analysis. 

Mantel path analysis 

Path analysis is a regression technique to explore causal connec- 
tions among relevant factors. We used Mantel path analysis to 
quantify direct and indirect influence of various biotic and abi- 
otic factors on spatial distribution and performance of rockcress 
population. Partial Mantel tests were conducted to evaluate the 
plausible hypotheses (generated by combining the information 
from correlograms and pairwise Mantel tests) on underlying 
ecological processes by looking at the relationship between two 
variables and keeping all other variables constant. For exam- 
ple, given the a priori hypothesis that A effects B and B effects 
C, the partial Mantel test can be used to test if C effects A in 
the absence of B, C~ A|B. If A and/or C displayed significant 
spatial autocorrelation then space was added as an additional par- 
tial, C~A|Space + B, in the partial Mantel test (See Table Al in 
Appendix for all the hypotheses). The direction of relationships 
was predetermined based on the common ecological knowledge. 
For instance, local relief influences soil moisture and not the other 
way. Similarly, the number of reproductive stalks influences the 
number of seeds and not vice versa. When the direction could not 
be determined using this approach then it was left as a simple cor- 
relation between two variables. The R package "ecodist" (Goslee 
and Urban, 2007) was used to calculate the Mantel correlogram 
and conduct Mantel path analysis. 

RESULTS 

WHAT GOVERNS THE SPATIAL DISTRIBUTION AND PERFORMANCE OF 
A ROCKCRESS POPULATION? 

Microhabitat preference 

Rockcress plants were not evenly distributed across the local 
landscape, occurring at highest densities at the lower end of 



the N-facing slope and at relatively high soil moisture levels 
(Figures 1A-F). All variables displayed significant spatial auto- 
correlation except soil moisture (Figures 2A-C). Elevation dis- 
played significant spatial autocorrelation (Figure 2A) and a neg- 
ative relationship with soil moisture (Mantel tests, P = 0.05, 
Figure 3A), which was also present after accounting for the spa- 
tial structure (partial Mantel tests, Figure 4). In general, rockcress 
size and reproduction showed a negative relationship with ele- 
vation, a positive relationship between soil moisture and no 
relationship with vapor pressure deficit. Thus, rockcress occurred 
more often in lower, moist areas and in these areas growth and 
reproduction were also higher. These patterns were reflected in 
correlograms showing significant (P < 0.05) spatial autocorrela- 
tion among rockcress individuals (Figure 2B) and Mantel tests 
showing significant negative correlation between elevation and 
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FIGURE 2 | Mantel correlograms showing significant spatial 
autocorrelation for (A) abiotic environmental factors (Elevation, and 
Vapor pressure deficit), (B) intra-specific performance and fitness 
measures [rosette diameter, stalk number, stalk diameter, fruit 
number, and line), and (C) inter-specific neighbors and herbivory 
{goldenrod density, dandelion density, rosette herbivory, and stalk 
herbivory). Filled (significant) and open (non-significant) circles represent 
statistical significance (P < 0.05) of Mantel correlation coefficients and 
each point marks the middle point of the respective lag distance bin. The 
range of autocorrelation for each variable is given next to the legend. Only 
significant (P < 0.05) correlograms are shown. 
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(B) intra-specific interactions, and (C) inter-specific interactions. Error 
bars mark the 90% confidence interval for correlation coefficients. 
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rockcress (rosette diameter, stalk number, and fruit number) as 
well as a positive relationship between soil moisture and rock- 
cress (stalk number) (Figure 3A). The relationship between soil 
moisture and stalk number was not significant after accounting 
for space (Figure 4). Vapor pressure deficit had no influence on 
rockcress size or fecundity. Dandelion density around rockcress 
was highest toward mid-elevations (Figure ID) and no influence 
of vapor pressure deficit was detected. The density of golden- 
rod plants around rockcress was positively correlated with high 
vapor pressure deficit (Figure 3 A) but this relationship was not 
significant after accounting for space. 

INTRA-SPECIFIC INTERACTIONS 

Rockcress growth (rosette diameter), reproduction (stalk num- 
ber, fruit number), and line were significantly spatially autocor- 
related (Figure 2B) with varying ranges from 6.4 m (line) to 
20.8 m (rosette diameter). Genetically similar individuals were 



found in close proximity of one another (spatially autocorre- 
lated, Figure 2B) and rosette diameter had positive influence on 
reproduction with and without space (Figures 3B, 4). Spatial dis- 
tribution patterns of rockcress had both genetic (Figure 2B) and 
environmental (Figures 3A, 4) components. 

INSECT HERBIVORY ON ROCKCRESS 

Herbivory on rockcress rosette and reproductive stalk leaves 
showed significant (P = 0.05) spatial autocorrelation [rosette 
herbivory (range = 12.6 m), stalk herbivory (range = 8.4 m)] 
(Figure 2C). Spatial cross-correlation (Mantel tests, Figure 3C) 
and cross-correlations after removing multiple correlation and 
space (Partial Mantel tests, Figure 4, Table Al in Appendix) 
showed positive relationship between rosette herbivory and stalk 
herbivory. Interspecific neighbor density (dandelion density) 
around rockcress was positively correlated with rosette herbivory 
(Figures 3C, 4), and herbivory on the leaves of reproductive stalks 
(stalk herbivory) was negatively correlated with fruit production 
(fruit number) (Figures 3C, 4), indicating negative impact of 
rosette herbivory on rockcress reproductive performance. 

INTER-SPECIFIC INTERACTIONS 

Mantel correlograms showed significant (P = 0.05) spatial auto- 
correlation in neighboring densities of goldenrod (range = 
12.6 m) and dandelion (range = 12.6 m) (Figure 2C). Partial 
Mantel tests did not indicate any relationship between rockcress 
and neighboring densities of goldenrod and dandelion. 

WHAT IS THE RELATIVE IMPORTANCE OF FACTORS GOVERNING THE 
SPATIAL DISTRIBUTION AND PERFORMANCE OF A ROCKCRESS 
POPULATION? 

The spatial distribution of rockcress individuals across the land- 
scape was greatly influenced by the spatial distance (significant 
spatial autocorrelation, Figure 2B) and local relief (Figure 3). 
Rockcress plants were generally crowded in N-facing slopes at 
lower elevations where there is shade by trees (Figure 1). Soil 
moisture and vapor pressure deficit did not show any signifi- 
cant influence on spatial distribution of rockcress plants. Range 
of autocorrelation for genotypic diversity (line) was the lowest 
among intra-specific traits that indicates greater intra-specific 
diversity of rockcress plants within short distances. Inter-specific 
interactions showed some indirect influence on rockcress distri- 
bution mediated through increased herbivory (Figure 4). Overall 
space and elevation and has the greatest influence on spatial 
distribution of rockcress individuals. 

Elevation, soil moisture, and vapor pressure deficit did not 
display any significant influence on rockcress performance. Intra- 
specific interactions and insect herbivory are the main drivers 
of rockcress performance. Additionally, indirect influence of 
inter-specific interactions on rockcress performance were evident 
(Figures 3, 4). 

DISCUSSION 
MICR0HABITAT PREFERENCE 

Spatial autocorrelation for plant performance and genetic varia- 
tion indicate spatial aggregation among related plants for general 
fitness, which is most likely a consequence of limited dispersal 
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FIGURE 4 | An ecological network of interactions between abiotic 
environmental factors (Elevation, Vapor pressure deficit. Soil 
moisture, and Nutrients), intra-specific performance and fitness 
measures [rosette diameter, stalk diameter, fruit number, stalk 
number, and line), and inter-specific neighbors and herbivory 
(goldenrod density, dandelion density, rosette herbivory, and stalk 
herbivory) in a rockcress population. Boxes represent observed 
variables and arrows represent significant partial correlations. Sign 



associated with each arrow describes the positive or negative influence. 
Thickness of arrows is related to the Partial Mantel correlation 
coefficient associated with each arrow and 90% confidence interval is 
given in the parentheses. Unidirectional arrows represent positive or 
negative influence and bidirectional arrows represent simple correlations. 
The dotted gray boxes display the variables not measured in this study 
and dotted gray arrows describe the proposed mechanisms to explain 
the relationships. 



and habitat heterogeneity. Facilitation among con-specific plants 
is a less likely explanation, given inter-plant spacing among rock- 
cress individuals was typically much greater than the reach of 
canopies or root systems; the nearest neighbors of the rockcress 
plants were often other species of plants. Genetic similarities 
decreased at shorter distances (<5m) than did performance or 
fitness measures (8-20 m) and landscape attributes (e.g., topog- 
raphy, ~16m), indicating that clustering of similar sized plants 
had both genetic and micro-environmental influences. To date, 
local adaptation in rockcress for defensive traits has been doc- 
umented on a broad geographic scale (Prasad et al., 2012). We 
propose that a greater diversity of functional genes could be 



present in the rockcress population at finer geographical scales 
due to adaptation to local microenvironments, limited dispersal, 
and a predominantly self-fertilizing breeding system (Song et al., 
2006). 

Elevations in the Black Hills (e.g., our study site at about 
1700 m) are much lower than elevations where this species 
(B. stricta) of rockcress is usually found in other mountain 
ranges (about 2500 m) (Mitchell-Olds, personal communica- 
tion). Probably as a consequence of lower and drier habitats 
in the Black Hills, we find rockcress on N-facing slopes and in 
shaded areas. Recent genotyping studies indicate that rockcress 
in the Black Hills probably originated from lower latitudes in the 
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southern Rockies (Lee and Mitchell-Olds, 2011) where rockcress 
may be adapted to drier climates. This study was focused on the 
spatial characteristics and only collected environmental data for 
1 day across the meadow. Further work is needed to assess the 
different spatial and temporal scales at which local adaptation, in 
particular to dry environments, may occur in this system. 

INSECT HERBIVORY ON ROCKCRESS 

Insect herbivores are known to reduce the fitness of plants in the 
wild (Marquis, 1992) and previous studies (Bloom et al., 2003; 
Prasad et al., 2012) on rockcress have reported a significant neg- 
ative correlation between rockcress fecundity and herbivory on 
rosette and reproductive stalk leaves. However, the direct and 
indirect environmental effects on susceptibility to herbivores are 
also likely to be important in this system. Based on significant spa- 
tial autocorrelation and strong positive correlation of flea beetle 
herbivory on rosette and reproductive stalk leaves of rockcress, we 
suggest that plants in close proximity to an infested plant are more 
prone to flea beetle attack than the plants that are further away. 
This spatial clustering of herbivory among plants could possibly 
be due to (1) direct and indirect microhabitat (indirect effects of 
the environment include effects on host plant resistance through 
stressful environments and competitors), (2) spatial association 
of plants within a line with similar susceptibilities and resis- 
tances, and (3) that flea beetles tend to attract one another. It 
should be noted that we did not detect a significant correlation 
between line and herbivory, despite the detection of significant 
genetic variation from common garden experiments among these 
lines in resistance and glucosinolate toxin production in previ- 
ous studies (Siemens et al., 2009). Thus, the micro-environmental 
component to spatial variation in herbivory appears to be more 
important than genetic variation. Additionally, the positive cor- 
relation between rosette herbivory and dandelion density indicates 
an environmental component to the susceptibility. Because rock- 
cress and dandelion plants did not often occur in close proximity 
to one another for any competitive interactions to occur between 
them, we suggest that increased densities of dandelion may be 
an indicator of a certain habitat quality that somehow affects 
rockcress susceptibility. For example, less optimal habitats in 
this system may stress plants and make them more suscepti- 
ble to herbivory. Molecular studies on Arabidopsis indicate that 
responses to drought stress may attenuate defense responses to 
disease and herbivores (Fujita et al., 2006). Additional stud- 
ies suggest that experimental drought stress in rockcress low- 
ers basal levels of glucosinolates (Alsdurf et al., 2013) and 
we have observed dramatic increases in herbivory on rock- 
cress plants during severe drought (David H. Siemens, personal 
observations). 

INTER-SPECIFIC INTERACTIONS 

Increased densities of dandelion and rockcress may nonetheless 
indicate increased competitive interactions because the meadow 
community studied was diverse (~41 plant species were recorded 
in 499 micro-communities [100 cm 2 ] around 499 rockcress indi- 
viduals in a common garden in this same meadow) and dense 
(~70 plants/ 100 cm 2 , n = 499) (Siemens and Haugen, 2013). In 
a common garden experiment in the study area set up across 



the local range boundary, Siemens and Haugen (2013) found 
that the decreased performance of rockcress was correlated with 
the change in community structure. Removal experiments in 
the field are needed to determine whether the stress that is 
correlated with community change is caused by competitive inter- 
actions (Mulder and Ruess, 1998; Van der Wal et al, 2000). 
Other species that thrive in similar elevation and soil moisture 
ranges to rockcress and goldenrod might be expected to be more 
competitive for space, nutrients, and moisture. However, rock- 
cress is likely to have weaker competitive interactions with other 
species, like dandelion, that occupy slightly different microhab- 
itats, i.e., relatively higher elevation and drier areas. Goldenrod 
did not have negative influence on rockcress while dandelion 
showed positive relationship with rockcress herbivory. Species 
like dandelion that occupy drier habitats may allocate more 
resources to roots and thus would be expected to be better 
below ground competitors. The root-shoot ratio of goldenrod 
[~0.4 (Johnson and Biondini, 2001)] makes it a potentially 
poorer competitor than dandelion [2.5 (Thomas and Bazzaz, 
1996)] and rockcress [1.5 (Haugen et al., 2008)]. Our results sup- 
port previous findings (Eskelinen, 2008) suggesting that positive 
and negative interactions may vary on target species. However, 
because of the significant spatial segregation, different habitat 
preferences, and overall low densities of dandelion and gold- 
enrod within the community it is more likely that the spatial 
associations and correlated effects on rockcress performance 
are caused by effects of the entire community, limited disper- 
sal and variation in microhabitats in space and time. Recent 
studies suggest the importance of intra- and interspecific geno- 
type interaction in structuring the plant community (Crutsinger 
et al., 2006; Lankau and Strauss, 2007; Genung et al., 2012) and 
results from this study indicate the spatial distribution has both 
genetic and environmental components. Overall, this study sug- 
gests that common microhabitat preference and limited dispersal 
are the main drivers for spatial structure of the rockcress popu- 
lation. However, intra-specific interactions and insect herbivory 
are the main drivers of rockcress performance in the meadow 
community. 

Abiotic and biotic interactions are interconnected strongly 
in space. The combination of complementary statistical tools 
(Bayesian kriging, the piecewise Mantel correlogram and Mantel 
path analysis) enhanced our understanding of the underlying eco- 
logical processes of complex spatial interactions and allowed us 
to dissect direct and indirect effects of biotic and abiotic factors 
on distribution and performance of a rockcress population. Our 
attention to spatial distributions and associations expands the 
possible causations underlying community assembly associations 
over non-spatial analyses and our results help generate hypothe- 
ses for future experimental studies in ecological and evolutionary 
genomics. 

ACKNOWLEDGMENTS 

This study was supported by NSF EAR-0405381, NSF/EPSCoR 
EPS-0091948, NIH/INBRE 2 P20 RR016479 grants. We thank 
Sarah Adelman and Riston Haugen for field assistance and Seth 
Wilberding for proving comments on the earlier version of this 
draft. 



www.frontiersin.org 



July 2014 | Volume 5 | Article 348 | 7 



Naithani et al. 



Spatial distribution and performance of Boechera striata 



REFERENCES 

Alsdurf, J. D., Ripley, T. J., Matzner, S. L.» and Siemens, D. H. (2013). Drought- 
induced trans-generational tradeoff between stress tolerance and defence: 
consequences for range limits? AoB Plants 5:plt038. doi: 10.1093/aobpla/ 
plt038 

Al-Shehbaz, I. A. (2003). Transfer of most North American species of Arabis to 
Boechera (Brassicaceae). Novon 13, 381-391. doi: 10.2307/3393366 

Arhonditsis, G. B., Stow, C. A., Steinberg, L, J., Kenney, M. A., Lathrop, R. 
C, McBride, S. J., et al. (2006). Exploring ecological patterns with struc- 
tural equation modeling and Bayesian analysis. Ecol. Model. 192, 385-409. doi: 
10.1016/j.ecolmodel.2005.07.028 

Becerra, J. X. (2007). The impact of herbivore -plant coevolution on plant 
community structure. Proc. Natl. Acad. Sci. U.S.A. 104, 7483-7488. doi: 
10.1073/pnas.0608253104 

Bloom, T. C, Baskin, J. M., and Baskin, C. C. (2003). Ecological life history of the 
facultative woodland biennial Arabis laevigata variety laevigata (Brassicaceae): 
effects of leaf litter cover, herbivory, and substrate-type on bolting and fecun- 
dity./. TorreyBot. Soc. 130, 16-22. doi: 10.2307/3557521 

Brooker, R. W., Maestre, E T., Callaway, R. M., Lortie, C. L., Cavieres, 
L. A., Kunstler, G., et al. (2008). Facilitation in plant communities: the 
past, the present, and the future. /. Ecol. 96, 18-34. doi: 10. Ill 1365- 
2745.2007.01295.x 

Callaway, R. M., and Walker, L. R. (1997). Competition and facilitation: a syn- 
thetic approach to interactions in plant communities. Ecology 78, 1958-1965. 
doi: 10.2307/2265936 

Crutsinger, G. M., Collins, M. D., Fordyce, J. A., Gompert, Z., Nice, C. C., and 
Sanders, N. J. (2006). Plant genotypic diversity predicts community struc- 
ture and governs an ecosystem process. Science 313, 966. doi: 10.1126/sci- 
ence. 1128326 

Diggle, P. J., Moyeed, R. A., andTawn, J. A. (1998). Model-based Geostatistics. Appl. 
Stat. 47, 299-350. 

Diggle, P. J., and Ribeiro, P. J. (2002). Bayesian inference in gaussian 
model-based geostatistics. Geogr. Environ. Model. 6, 129-146. doi: 
10.1080/1361593022000029467 

Eskelinen, A. (2008). Herbivore and neighbour effects on tundra plants depend on 
species identity, nutrient availability and local environmental conditions. /. Ecol. 
96, 155-165. doi: 10.1111/j.l365-2745.2007.01322.x 

Fujita, M., Fujita, Y., Noutoshi, Y., Takahashi, F., Narusaka, Y., Yamaguchi- 
Shinozaki, K., et al. (2006). Crosstalk between abiotic and biotic stress 
responses: a current view from the points of convergence in the stress sig- 
naling networks. Curr. Opin. Plant Biol. 9, 436-442. doi: 10.1016/j.pbi.2006. 
05.014 

Genung, M. A., Bailey, J. K., and Schweitzer, J. A. (2012). Welcome to the neigh- 
bourhood: interspecific genotype by genotype interactions in Solidago influence 
above- and belowground biomass and associated communities. Ecol. Lett 15, 
65-73. doi: 10. 1 1 1 1/j. 1461-0248.201 1.01710.x 

Goslee, S. (2010). Correlation analysis of dissimilarity matrices. Plant Ecol 206, 
279-286. doi: 10.1007/sll258-009-9641-0 

Goslee, S. C, Niering, W. A., Urban, D. L., and Christensen, N. L. (2005). Influence 
of environment, history and vegetative interactions on stand dynamics in 
a Connecticut forest. /. Torrey Bot. Soc. 132, 471-482. doi: 10.3159/1095- 
5674(2005)132[471:IOEHAV]2.0.CO;2 

Goslee, S. C, and Urban, D. L. (2007). The ecodist package for dissimilarity-based 
analysis of ecological data. /. Stat. Softw. 22, 1-19. 

Haugen, R., Steffes, L., Wolf, J., Brown, P., Matzner, S., and Siemens, D. H. 
(2008). Evolution of drought tolerance and defense: dependence of tradeoffs 
on mechanism, environment and defense switching. Oikos 117, 231-244. doi: 
10.1111/j.2007.0030- 1299.16111.x 

Holzapfel, C., and Mahall, B. E. (1999). Bidirectional facilitation and interference 
between shrubs and annuals in the Mojave desert. Ecology 80, 1747-1761. doi: 
10.2307/176564 

Jeltsch, E, Moloney, K., and Milton, S. J. (1999). Detecting process from snap- 
shot pattern: lessons from tree spacing in the Southern Kalahari. Oikos 85, 
451-466. 

Johnson, H. A., and Biondini, M. E. (2001). Root morphological plasticity and 
nitrogen uptake of 59 plant species from the Great Plains grasslands, U.S.A. 
Basic Appl. Ecol. 2, 127-143. doi: 10.1078/1439-1791-00044 

Knight, C. A., Vogel, H., Kroymann, J., Shumate, A., Witsenboer, H., and 
Mitchell- Olds, T. (2006). Expression profiling and local adaptation of Boechera 



holboellii populations for water use efficiency across a naturally occurring 
water stress gradient. Mol. Ecol. 15, 1229-1237. doi: 10.1111/j.l365-294X.2006. 
02818.x 

Lankau, R. A., and Strauss, S. Y. (2007). Mutual feedbacks maintain both genetic 

and species diversity in a plant community. Science 317, 1561-1563. doi: 

10.1126/science.ll47455 
Leduc, A., Drapeau, P., Bergeron, Y, and Legendre, P. (1992). Study of spatial com- 
ponents of forest cover using partial Mantel tests and path analysis. /. Veg. Sci. 

3, 69-78. doi: 10.2307/3236000 
Lee, C.-R., and Mitchell-Olds, T. (2011). Quantifying effects of environmental 

and geographical factors on patterns of genetic differentiation. Mol. Ecol. 20, 

4631-4642. doi: 10.1111/j.l365-294X.2011.05310.x 
Levin, S. A. (1992). The problem of pattern and scale in ecology: the Robert H. 

MacArthur award lecture. Ecology 73, 1943-1967. 
Loranty, M. M., Mackay, D. S., Ewers, B. E., Adelman, J. D., and Kruger, E. L. 

(2008). Environmental drivers of spatial variation in whole-tree transpiration 

in an aspen-dominated upland-to-wetland forest gradient. Water Resour. Res. 

44:W02441. doi: 10.1029/2007WR006272 
Mantel, N. (1967). The detection of disease clustering and a generalized regression 

approach. Cancer Res. 27, 209-220. 
Marquis, R. J. (1992). "The selective impact of herbivores," in Plant Resistance 

to Herbivores and Pathogens: Ecology, Evolution, and Genetics, eds R. 

S. Fritz and E. L. Simms (Chicago, IL: University of Chicago Press), 

301-325. 

Mclntire, E. J. B., and Fajardo, A. (2009). Beyond description: the active and 
effective way to infer processes from spatial patterns. Ecology 90, 46-56. doi: 
10.1890/07-2096.1 

Mitchell-Olds, T. (2001). Arabidopsis thaliana and its wild relatives: a model system 
for ecology and evolution. Trends Ecol Evol. 16, 693-700. doi: 10.1016/S0169- 
5347(01)02291-1 

Mulder, C. P. H., and Ruess, R. W. (1998). Effects of herbivory on arrowgrass: inter- 
actions between geese, neighboring plants, and abiotic factors. Ecol. Monogr. 68, 
275-293. doi: 10.2307/2657204 

Prasad, K. V. S. K., Song, B. H., Olson-Manning, C., Anderson, J. T., Lee, C. R., 
Schranz, M. E., et al. (2012). A gain-of-function polymorphism controlling 
complex traits and fitness in nature. Science 337, 1081-1084. doi: 10.1126/sci- 
ence.1221636 

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population 
structure using multilocus genotype data. Genetics 155, 945-959. 

Pugnaire, F. I., and Luque, M. T. (2001). Changes in plant interactions along 
a gradient of environmental stress. Oikos 93, 42-49. doi: 10.1034/j.l600- 
0706.2001.930104.x 

R Development Core Team (2012). R: A Language and Environment for Statistical 
Computing. Vienna: R Foundation for Statistical Computing. Available online 
at: http://www.R-project.org 

Ribeiro, P. J. Jr., and Diggle, P. J. (2001). geoR: A package for geostatistical analysis. 
RNews 1, 14-18. 

Rushworth, C. A., Song, B.-H., Lee, C.-R., and Mitchell-Olds, T. (2011). Boechera, 
a model system for ecological genomics. Mol. Ecol 20, 4843-4857. doi: 10.1 1 1 1/ 
j.l365-294X.2011.05340.x 

Siemens, D. H., and Haugen, R. (2013). Plant chemical defense allocation con- 
strains evolution of tolerance to community change across a range boundary. 
Ecol Evol. 3, 4339^347. doi: 10.1002/ece3.657 

Siemens, D. H., Haugen, R., Matzner, S., and Vanasma, N. (2009). Plant chemical 
defence allocation constrains evolution of local range. Mol Ecol 18, 4974-4983. 
doi: 10. 1 1 1 l/j.l365-294X.2009.04389.x 

Smithwick, E. A. H., Naithani, K. J., Balser, T. C, Romme, W. H., and Turner, 
M. G. (2012). Post-Fire spatial patterns of soil nitrogen mineralization 
and microbial abundance. PLoS ONE 7:e50597. doi: 10.1371/journal.pone. 
0050597 

Song, B.-H., Clauss, M. J., Pepper, A., and Mitchell -Olds, T. (2006). Geographic 
patterns of micros at ellite variation in Boechera striata, a close relative 
of Arabidopsis. Mol. Ecol. 15, 357-369. doi: 10.1111/j.l365-294X.2005. 
02817.x 

Song, B.-H., and Mitchell -Olds, T. (2011). Evolutionary and ecological 
genomics of non-model plants. /. Syst. Evol. 49, 17-24. doi: 10. 1 1 1 1/j. 1759- 
6831.2010.00111.x 

Thomas, S. C, and Bazzaz, F. A. (1996). Elevated C02 and leaf shape: are 
dandelions getting toothier? Am. J. Bot. 83, 106-111. doi: 10.2307/2445962 



Frontiers in Plant Science | Functional Plant Ecology July 2014 | Volume 5 | Article 348 | 8 



Naithani et al. 



Spatial distribution and performance of Boechera stricta 



Urban, D., Goslee, S., Pierce, K., and Lookingbill, T. (2002). Extending community 
ecology to landscapes. Ecoscience 9, 200-212. 

Van der Wal, R., Egas, M., van der Veen, A., and Bakker, J. (2000). Effects of resource 
competition and herbivory on plant performance along a natural productivity 
gradient./. Ecol. 88, 317-330. doi: 10.1046/j.l365-2745.2000.00450.x 

Wagner, H. H. (2003). Spatial covariance in plant communities: integrating 
ordination, Geostatistics, and variance testing. Ecology 84, 1045-1057. doi: 
10.1890/0012-9658(2003)084%5B1045:SCIPCI%5D2.0.CO;2 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 



Received: 07 March 2014; accepted: 28 June 2014; published online: 22 July 2014. 
Citation: Naithani KJ, Ewers BE, Adelman JD and Siemens DH (2014) Abiotic and 
biotic controls on local spatial distribution and performance of Boechera stricta. Front. 
Plant Sci. 5:348. doi: 10.3389/fpls.2014.00348 

This article was submitted to Functional Plant Ecology, a section of the journal 
Frontiers in Plant Science. 

Copyright © 2014 Naithani, Ewers, Adelman and Siemens. This is an open-access 
article distributed under the terms of the Creative Commons Attribution License 
( CC BY). The use, distribution or reproduction in other forums is permitted, provided 
the original author(s) or licensor are credited and that the original publication in this 
journal is cited, in accordance with accepted academic practice. No use, distribution or 
reproduction is permitted which does not comply with these terms. 



www.frontiersin.org 



July 2014 | Volume 5 | Article 348 | 9 



Naithani et al. 



Spatial distribution and performance of Boechera striata 



APPENDIX 



Table A1 | Hypotheses derived from correlograms and simple Mantel tests, and tested by using partial Mantel tests. 
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Bold p-value indicates significant correlation. 
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FIGURE Al | Semivariogams (A-C), prior and posterior distribution of range parameter (range = 34>) (D-F), and cross-validation graphs (G-l) 
associated with elevation, soil moisture, and vapor pressure deficit. 
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